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Abstract 

We present a scheme to prepare a quantum state in a ion trap with probability 
approaching to one by means of ion trap quantum computing and Grover's 
quantum search algorithm acting on trapped ions. 
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Quantum computing [1-5] is a field at the junction of modern physics and theoretical 
computer science. Practical experiments involving a few quantum bits have been successfully 
performed, and much progress has been achieved in quantum information theory, quantum 
error correction and fault tolerant. Although we still far from having desktop quantum 
computers in our offices, the quantum computational paradigm could soon be more than 
mere theoretical exercise. 

The discovery of Shor's quantum algorithm [6,7] for factoring whole numbers into prime 
integer factors was a major milestone in the history of quantum computing. Another signif- 
icant quantum algorithm is Grover's quantum search algorithm [8]. Grover's algorithm does 
not solve NP-complete problems in polynomial-time, but the wide range of its applications 
compensates for this. The key of performing quantum computations is the implementation 
of a quantum computer, because a quantum computer must meet the requirements: the iso- 
lation from the environment, the controllable and coherent interaction between the qubits, 
and the high-efficiency interrogation of individual qubits. Cirac and Zoller have proposed an 
attractive scheme for realizing a quantum computer, which is easily scalable to an arbitrary 
number of qubits. Their scheme is based on a collection of trapped atomic ions, where each 
qubit is comprised of a pair of the ions' internal states. The trapped-ion system features 
long coherent times of the internal qubit states, while the strongly-coupled ion motion al- 
lows quantum information to be transferred between different ions by using a particular 
quantized mode of the ions' collective motion. The crucial step is that the "quantum data 
bus" must be initialized in its ground state. The basics of this scheme have been demon- 
strated experimentally in a fundamental logic gate operating between a motional mode of 
a single trapped ion and two of the ion's internal states[9]. Recently, B.E.King et al. have 
successfully realized cooling the collective modes of motion of two trapped ions into their 
ground state[10]. This research pushed the quantum computation by using trapped ions to 
a new level. In this rapid communication, we show that a variety of multi-phonon coherent 
states in a ion trap can be efficiently prepared by combining ion trap quantum computing 
and quantum search algorithm. 
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We consider a register which consists of m two-level ions forming a string confined in a 
linear rf trap. The electronic states of ions are denoted by an integer k, 

\k) e = \s m ) <S> \s m -i) <g> • • • <S> \si) (1) 

where m is the number of the qubits in the register, Sj = 0, 1 and 

k = s m x 2 m ~ l + s m _! x 2 m ~ 2 + • • • + Sl x 2° (2) 

The interaction between the center-of-mass motion and the internal electronic state of 
each ion was mediated by a detuned laser pulse. For interaction times much greater than 
the vibrational period of the trap, the effective Hamiltonian for the jth ion [11] is 

H\ j) = M + a X (a^ + \) (3) 

where is the population inversion for jth ion, x = f] 2 ^i 2 /{mA) with r\ being the Lamb- 
Dicke parameter, Q being the Rabi frequency for the transitions between the two internal 
states of the ions, and A being the detuning between the exciting laser pulses and electronic 
transition. We assume that these parameters are the same for all ions. We choose the 
durations tj of the standing wave pulse to satisfy x r j — 2%/2 m , then the Hamiltonians 
generate a unitary transformation 

~ 2nia + aj 

u = ex P ( ^r~) ( 4 ) 

where the electronic operator 7 is defined by 

m 1 JV— 1 

7 = E(^ j) + 9) 2 '" 1 = E k \ k )M (5) 

3=1 Z k=0 

where N = 2 m is the size of the register space. 

Let $ is another electronic operator which is canonically conjugate to 7, the relation 
between the eigenstate \p) e of <E> and the eigenstate | k) e of 7 is connected by the following 
Fourier transformation 

Y N-l 

\P)e = -7= E ex P {-^ikpN- v )\k) e (6) 



where < p < N - 1. 

We assume that all ions are in the ground state |0) e and the vibrational state is arbitrary. 
Then the initial state of the system is 

oo 

\An) = Y; C n\n)vib® |0) e (7) 
n=0 

Now we perform three steps on above equation: Firstly we apply a sequence of 7r/2- 
pulses on the ions to produce a uniform superposition over all possible electronic energy 
eigenstates, which is also the state |0) e in the Fourier transform basis. Secondly, the unitary 
interaction in Eq.(4) is implemented to couple the vibrational and electronic states (note 
that U does not affect the vibrational state, but only displaces the eigenstates of the operator 
$). Finally, an inverse Fourier transform is run on the electronic register. After these three 
steps, the system approaches the following state[12] 

oo JV-1 

|Vw) = EE C k+nN\k + nN) vib <g> \k) e (8) 

n=0 k=0 

If the vibrational state of the ions is a coherent state 

\a\ 2 00 a n 

\ip)vib = \<x)vib = exp ( — ) Y ~^\n)vib, (9) 

I n=0 Vn! 

then Eq.(8) becomes 

N-1 

\i>out) = E W,N,k) vib ®\k) e (10) 

fc=0 

where 

1 00 a k+nN 

\a,N,k) vib = exp(--|a| 2 ) V =\k + nN) vib (11) 

2 „=o y(HniV)! 

Eq.(10) is an entangled state, the dimension of the entangled space is N. If we repeatly 
perform a series of Grover algorithm on Eq.(10), then the probability of finding one certain 
state \a, N,ko)vib ® \ko)e will becomes very large, and all the others probability become 
approximately zero. In order to do these, we rewrite Eq.(10) as 

N-1 

\ipout) = \a,N,k )vib® \k )e+ Yl \oc,N,k) vib ®\k) e (12) 



Then we perform the Grover's algorithm [13] on Eq.(12) by two steps. Firstly, we rotate 
the marked state \a, N, k ) vib <S> \ko) e by a phase of n to make Eq.(12) becomes the following 
form 

\ipout) = -\a,N,ko) vib ® \k ) e + J2 \u,N,k) vib ®\k) e (13) 

k=0,k^k 

Secondly, we continuously apply the inversion operation about average amplitude to 
Eq.(13) by the following unitary matrix 



Dij 



f -1 if < = j 



(14) 



After performing the above Grover's algorithm one time, we get 



JV-l 

\ip<mt)i = ai|a, N,k ) vib <g> \k ) e + b 1 \a,N,k) vib <S> \k) e (15) 

k=0,k^k 

where a x = VN sin 30, h = y/(N - 1)/Ncos30 with sin# = 1/^N. 

If we repeat the above Grover's algorithm j times on Eq. (12), we can obtain that 

\ipout)j = aj\a,N, k ) vib ® \k ) e + bj\a, N, k) vib <g> \k) e (16) 

k=0,k^ko 



where a 6 = v / A r sin[(2j + 1)0] and bj = ^J(N - l)/JVcos[(2j + 1)9]. 

From Eq.(16) we can see that if we repeat Grover's algorithm T = (n — 29) /AO times, 
then we have br = 0. If we now make a measurement on quantum states of the register, we 
can gain the state \ko) e with probability approaching to one. After the measurement, the 
vibrational state of ions becomes 

1 00 a ko+nN 

\a, N, k ) vib = exp ( — — |ct | 2 ) V . =\k + nN) vib (17) 

2 n=o ^J(k + nN)\ 

The above equation is just the fc -th eigenstate of operator a N . If we just want to get 
a Schrodinger-cat state, we need only one ion (m = 1, N = 2). In fact, when t ko — 0, 
we can get the even Schrodinger-cat state \a) + | — a), when ko = 1, we can get the odd 
Schrodinger-cat state \a) — \ — a). If we want to get four-phonon coherent states, we need 



only two ions, i.e., m — 2, and N — 4. By selecting /c = 0, 1, 2, 3, we can get the following 
four-phonon coherent states[14]: 

— \a) + \ — a) + \ia) + | — ia) (ko = 0) (18) 
1^)2 — \a) — \ — a) — i\ia) + i\ — ia) (k = 1) (19) 
|a) 3 — \a) + \ — a) — \ia) — \ — ia) (k = 2) (20) 

|a)4 = \a) — \ — a) + i\ia) — i\ — ia) (ko = 3) (21) 

It is clear that by means of this method, we can get all the eigenstates of operator a 2 ™ 
(m being the number of ions). 

In summary we provide a method which can produces multi-phonon coherent states 
with a completely successful probability. If we do not use Grover's algorithm, but measure 
immediately the electronic state in Eq.(12), then we obtain the state \k ) e with only a very 
less probability when iV is very large. We must note that Grover's algorithm can only repeat 
integer times, but T = (ir — 26)/ AO is not always an integer. Therefore we can only repeat 
Grover's algorithm by an integer nearest to T. Thus bx is not completely zero, we get the 
electronic state \k ) e with probability only approaching to one but not exact one. 

This work was supported by NSF of China, the Excellent Young- Teacher Foundation of 
the Educational Commission of China, ECF, STF and NSF of Hunan Province. . 
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